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Abstract 

The total neutron-Nucleus cross section has been calculated within an approach which 
takes into account nucleon-nucleon correlations, Glauber multiple scattering and in- 
elastic shadowing corrections. Nuclear targets ranging from 4 He to 208 Pb and neutron 
incident momentum ranging from 3 to 300 GeV/c, have been considered. Correlations 
have been introduced by two different approaches leading to the same results. The 
commonly used approximation, consisting in treating nuclear effects only by a product 
of one-body densities, is carefully analyzed and it is shown that the effects of realistic 
correlations resulting from modern nucleon-nucleon interactions and realistic corre- 
lations resulting from realistic nucleon-nucleon interactions and microscopic ground 
state calculation of nuclear properties cannot be disregarded. 
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The total neutron-Nucleus cross section at high energies has been the object of many 
calculations for its dependence is very sensitive to various effects, such as Glauber elastic 
PQ) and Gribov inelastic [2J) diffractive shadowing, which are relevant for the interpretation 
of color transparency phenomena [2111] and relativistic heavy ion processes [5]. The major 
mechanism governing the total cross section is Glauber inelastic shadowing, but a quantita- 
tive explanation of the experimental data has been achieved in the past only by considering 
also the effects of inelastic shadowing [6l [TJ. All calculations so far performed were based 
upon the so called one body density approximation, in which all terms but the first one of 
the correct expansion of the square of the nuclear wave function in terms of density matrices 
[5] are disregarded, which amounts to neglect all kinds of nucleon nucleon correlations. The 
necessity to investigate the effects of correlations on the total cross section a l °\ was pointed 
out by several authors pi [7J. It is precisely the aim of this work to present the results 
of calculations of the total neutron-Nucleus cross section within an approach based upon 
realistic many-body correlated wave functions (9] obtained with realistic nucleon-nucleon 
interactions [10J, Glauber multiple scattering and Gribov inelastic shadowing. 



1 Basic formalism 

Considering both Glauber (G) elastic scattering and Gribov inelastic shadowing (IS) , the 
total cross section on a nucleus A can be written as follows 

4-7T 

of = ai + a{ s = T Im [*£(0) + F«(0)] (1) 

where and a A denote respectively the Glauber and inelastic shadowing contributions, 
and -Foo(O) and -^(0) the corresponding forward elastic scattering amplitudes related to 
the full nuclear profile r^ 5 '' as follows 

^ <,s, (o) = |/*„rS' s, (M (2) 

The Glauber nuclear profile describing the elastic scattering of the neutron has the usual 
form 

A 

rfo(&n) = 1 " n^lt 1 - r N(Pn-»jWo), (3) 

where ipo = T/>o( r i> r 2, i*3> •••I'a), with r,- = (sj,Zj), is the ground state wave function of 
the target nucleus, b n the impact parameter of the neutron moving along the z-axis, and 
-T/v(6 n ) the nucleon-nucleon elastic profile function. The Inelastic Shadowing Profile should 
describe the diffractive dissociation of the neutron via the process n + iV — > X + N and 
its de-excitation to the ground state by the process X + iV — >• n + N, as well as the elastic 
scattering of X off the target nucleons. The three processes are described by the inelastic 
profiles T NX and T X n, and by the elastic profile Fx, respectively. In our approach, as in 
Ref. P|, we will consider only two non-diagonal transitions, i.e. n + iV — ► X + iV and 
X + iV — * n + N, and the elastic scattering of X. The corresponding diagrams are shown 
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in Fig. [H Within such an approximation, the Inelastic Shadowing profile can be written in 
the following form [3]: 



X I i<j 

A 

x Yl i 1 - T x{K - b h )\ Q(z k - Zi ) e(zj - z k ) x 

x [J [1 - r w (^ - 6,)] 6(zi - «,) e( z , - z 3 .) w J (4) 

where 

Qx = K ~ k x (5) 

is the longitudinal momentum transfer. The basic nuclear ingredient appearing in Eqs. ([3]) 
and (j3j) is the square of the nuclear wave function ip , which can be written in terms of 
density matrices as follows [8]: 

A 

\ip (r 1 ,...,r A )\ 2 = Yl Pi( r j) + ^2 A ( r u r j) Yl Pi( r k) + 

j=l i<j kj^i,j 

+ ^A(r i) r i )A(r fe) r 1 )[[ ft (r m ) + ...; (6) 

(i<j)yt(k<l) m=£i,j,k,l 

in which pi{rj) is the one-body density matrix 

r A 

Pi(n) = / \Mri,-,r A )\ 2 s(j2ri) (7) 

J j=2 

and the two-body contraction A is defined as follows: 

A(t*i,7* 2 ) = p 2 {ri,r 2 ) - pi{r x )p 1 {r 2 ) , (8) 
where P2{fi, 1*2) is the two-body density matrix 

r A 

Yldr, \Mri,r 2 ...,r A )\ 2 6(j2(r J )) ■ (9) 

J j=3 

The one- and two-body density matrices appearing in Eq. ([6]) are normalized according to 

J drpx(r) = J dr 1 dr 2 p 2 (r 1 ,r 2 ) = 1 (10) 

and satisfy the following sequential conditions: 

J drjp x {r h rj) p x (rj,r k ) = pi{r h r k ) (11) 

3 



which leads to 



dr 2 p 2 (r 1 ,r 2 ) = Pi{r x ) 



Jdr 1 A{r 1 ,r 2 ) = J dr 2 A(n, r 2 ) = 0; 



(12) 



(13) 



In Eq. ([6]) only unlinked contractions have to be considered, and the higher order terms 
include unlinked products of 3, 4, etc. two-body contractions, unlinked products of three- 
body contractions, describing three-nucleon correlations, and so on. When all terms up to 
A-body correlations are written down explicitly, an identity is obtained. 

The common approximation in Glauber type calculations consists in disregarding all 
terms of Eq. ([6]) but the first one. In this case the very well known expression for the total 
Glauber profile is given by 

A 



r G (&n) 



1 



1 



dr 1 p 1 (ri) T(b n - si) 



(14) 



By taking into account two-body correlations, i.e. all unlinked products of two-body con- 
tractions in Eq. (Q, one obtains [TT| [T2] 

A 



1 - 



X 



4-^ 



E 

771=0 



drxp^rx) T(b n - 
,i! |1 Jdr l dr 2 A(r u r 2 )T(b n -s 1 )T(b n -s 2 ) 



{A-2m)\m\ 2 



[l- J dr 1 p 1 r 1 )T(b n - s 1 



(15) 



which in the optical limit (A >> 1) becomes 

-A J dr x p x {r x ) T(b n - s x ) + ^ J dr x dr 2 A(n, r 2 ) Y{b n - 8l ) T(b n - s 2 ) 



(16) 

As for the inelastic shadowing contribution (J4J), it can be reduced to an expression depending 
upon the total nucleon and diffractive cross sections crff and <t*°* respectively. 

47T 

T 



a 



IS 



ImT^{0) = 2 J db n Til{b n ) = 

db n dz 1 dz 2 p(b n ,z 1 )p(b n ,z 2 ) / dM 2 x 

J d q TdMx 

. ' 0( Zl _ Z2 ) e i(Plab ~ Pm)(z 2 - Z X ) . 

• e- A / dz' p(b n , z') £ ■ e~ A / dz' p(b n , z') ^ ■ e~ A / dz' p(b n , z') ^ + 

Jz 2 Z J-oo Z JZi Z 

+ e ( Z2 _ Zi y(piab-Pm)(zi-Z 2 ). 

fZ\ l>Z\ P+CO 

■ e~ A / dz' p(b n , z>) £ ■ e~ A / dz' p(b n , z') ^ ■ e~ A / dz' p(b n , z') ^ + 

Jz 2 1 J-oo 1 J Z 2 1 . 

d Zl dz 2 e 1 ^ ~ **) A(b n , Zl ,b n ,z 2 ). (17) 



q T =0 

-+oo 
'Zi 



Within the approximation aff = a* oi and disregarding correlations (A(b n , z%, b n , z%) = 0) 
the well-known Karmanov-Kondratyuk |T3] expression is obtained 



i 8 = -An A 2 / db 



where depends upon M, 



d 2 a 



d 2 q T dM\ 



— — T(b) 

dM 2 x e 2 { >\F( qL ,z)\ 2 , (U 

q T =0 



In our calculations we have used both Eqs. ( |T6l) and (fTBl) . We have checked that the 
optical limit for A = 16 is valid within 1%, whereas correlations produce very tiny effects 
on a A s . The ingredients of our calculations were as follows: 

1. The density matrices have been obtained by a linked cluster expansion for the one- 
and two-body density operators expectation value, evaluated over a fully-correlated 
wave function [9] obtained variationally with the Argonne V8' interaction [10]. The 
one-body density has been obtained by integrating the two-body density. Let us stress 
that, unlike previous calculations, our two-body contractions (Eq. (jBJ)) exactly satisfy 
the condition given by Eq. (TT51) : 



2. the Glauber profile function is of the usual form 



r ( b ») = (1 - ia) e" 6 ^ (19) 

4:71 Oq 



with the energy- dependent parameters taken from [14J; 
3. the parameters for the inelastic shadowing were taken from [B]. 

The results of calculations for A He, 12 C , 16 and 208 Pb are presented in Fig. [2j The left 
panel shows the results obtained without correlations, whereas the effects of correlations 
are presented in the right panel. The results presented in Fig. [2] deserve the following 
comments: 

1. it can be seen that correlations increase the total cross section by about 10%, i.e. they 
decrease the nuclear transparency, worsening the agreement with the experimental 
data when only Glauber shadowing is considered; the inclusion of inelastic shadowing 
brings back theoretical calculations in good agreement with experimental data; 

2. in the case of A He we have calculated the cross section to all orders of correlations 
using the exact wave function of Ref . [15J ; it turns out that three- and four-nucleon 
correlations produce negligible effects on the total cross section; 

3. the effect of correlations is of the same order as the one from inelastic shadowing. 

It should be pointed out that the contribution to the optical phase shift (the second 
term in the exponent of Eq. ffl6|) ) is always negative; the black disk limit of our approach 
is satisfied. 
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In Fig. [3] the difference between the two-body density calculated within the mean field, 
the cluster expansion and the following approximation for the two-body density 



0j(ri,r a ) = pi{ri) pi{r 2 ) g{\r 1 - r 2 \) (20) 

frequently used in case of complex nuclei (see e.g. Ref. [H]), is exhibited. The curves 
represent the quantity 

p 2 (r) = j ' dRp 2 (r 1 = R+\r , r 2 = R-^rj . (21) 

It can be seen that the two-body density of Eq. (T2T]) represents a poor representation of the 
realistic one. We have checked to what extent the sequential relation (fT3|) is violated by the 
approximate two-body density matrix. The amount of violation of the integral in Eq. (TT3]) 
can be checked by calculating the quantity 

S(r 2 ) = J d ri A(r u r 2 ); (22) 

which is shown in Fig. H] for various nuclei. Fig. [5] shows the effect of violation of the 
sequential relation on the total cross section. 



2 A cluster expansion approach to the total cross sec- 
tion 

We have developed a cluster-expansion [9] formulation for a l °\ based upon the one-body dis- 
torted density matrix of Ref. [16] , obtained taking into account two-nucleon correlations at 
first order of the cluster expansion and Glauber rescatterings at all orders. The zeroth-order 
approximation (i.e. with no correlation effects) is the same as Eq. (TH1) ; correlations can 
be included with the first term of the wave function expansion of Eq. ([HD , by replacing the 
one-body densities appearing in such a term with the distorted one-body density proposed 
in Ref. [16], in such a way one obtains contributions representing the interaction of the 
incident nucleon with the particles involved in each of the diagrams contributing to the dis- 
torted density (namely, the shell model, hole and spectator diagrams). The final expression 
for the total cross section reads as follows: 

< 0t A = + + Aaf , (23) 
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in which the shell model (SM), hole (if) and spectator (S) contributions are as follows: 



a 



tot 
SM 



dbr, 



db r 



1 - 



1 

4 
A 



4 
A 



dr 1 p (r 1 )T(b n ,b 1 ] 



A tot 
a S 



dr l p (r l )T(b n ,b 1 
X ~A / dVldr2 H 12 Po{r\) Po{r 2) ~ Hf 2 \ Po (r 1 ,r 2 
-2 J db n (l - i J dr lPo {r x )Y{b ni b x )\ x 

y dr 1 dr 2 dr 3 p (r 1 ,r 2 )(AH^p (r 2 ,r 1 ) p a (3) + 

- ff| Po (r 2 , r 3 ) p (r 3 , n)) G^) G 2 (6 n ) r^) . 



[Gi(6 n )G 2 (6 n ) - 1] 



(24) 



Note that in the final calculations the shell model term has been exponentiated as in Eq. 
(USD- 

The results of calculations for 16 obtained with Eq. f[2~4"|) are compared in Fig. \5[ with 
the results predicted, by Eqs. (|SJ) and (l20j) . respectively. It can be seen that the expansion 
based on the distorted density of Ref. [16] is in perfect agreement with the one of the 
approach based on the expansion of the wave function, Eq. (flEl) . despite the different class 
of diagrams appearing in each contribution. 

It should be stressed that the expansion used for the distorted one-body density can 
be used to calculate, as in Ref. [IB], the total transparency in A(e,e'p)X experiments 
and distorted momentum distributions, as in Refs. ( [T71 [18]) with the full correlated wave 
function for complex nuclei. Calculations of nuclear and color transparencies in (e, e'p) and 
(p, 2p) are in progress and will be reported elsewhere. 



3 Summary and conclusions 

We have developed a method which can be used to calculate scattering processes at medium 
and high energies within a realistic and parameter-free description of nuclear structure. Our 
calculations followed the following strategy: 

i) the values of the parameters pertaining to the correlation functions and the mean field 
wave functions, have been obtained from the calculation of the ground-state energy, 
radius and density of the nucleus using realistic nucleon-nucleon interactions; 

ii) using these parameters we have calculated the total neutron-Nucleus cross section 
taking rigorously into account two-nucleon correlations within the expansion of the 
exact wave function [61 We have also adopted a cluster expansion procedure [16] 
obtaining essentially the same results. This gives us confidence that the treatment of 
correlations is model independent to a large extent. 
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The method we have developed appears to be a very effective, transparent and parameter- 
free one and the main results we have obtained are: 

i) the effects generated by two-nucleon correlations (i.e. by those parts of the wave 
function expansion ([6]) which contain two-body contractions), are of the same order as 
Gribov inelastic shadowing; this, in our opinion, points to the necessity of an analysis 
of the accuracy of the common approximation used in medium-high energy hadronic 
scattering processes consisting in disregarding all terms of the expansion ([6]) but the 
first one; 

ii) correlations due to three and higher order contractions appear to produce only negli- 
gible effects on the total cross section. 
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Figure 1: Typical diagrams describing elastic neutron- nucleus scattering at high energies. 
a) Glauber multiple scattering; b) and c) inelastic shadowing. 
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Figure 2: The neutron-nucleus total cross section, for 4 ife, 12 C, 16 and 208 Pb. Left panel: 
the result without the inclusion of correlations; dotted curves: one-body (first) term of Eq. 
(JTBJ) ; dot-dashed curves: one-body term plus inelastic shadowing effects (Eq. [TH]) . Right 
panel: results with the inclusion of correlation. Dotted curves: one-body (first) term of Eq. 
ffT6l) ; dashed curves: one-body term plus two-nucleon correlations of Eq. f|T6|) ; solid curves: 
one-body term plus two-nucleon correlations of Eq. ( |T6l) plus inelastic shadowing effects 
(Eq. [TBI) . Experimental data from [TT1 [T9] 
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Figure 3: The two-body density matrix of 16 within various approximations: the mean 
field approximation (dots), the approximation of Ref. [H] (dashes), the realistic calculation 
of the present work (full). 
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Figure 4: The quantity (J22l) . calculated within the approximation P2{fi-,i'2) = 
Pi( r i) Pi( r 2) 9(^12), for various nuclei. 
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Figure 5: Results of calculations for a to t using the pr>(r,r') of Refs. jT6], [TTj (red, dashed 
curve) as compared to the result corresponding to the wave function expansion of Ref. [8] 
(black, solid curve); the blue, dash-dotted curve is obtained with the correlation model of 
Ref. ITT] the black, dotted curve is the usual mean- field result. 
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